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Abstract 

The inert doublet model, a minimal extension of the Standard Model by a second higgs 
doublet with no direct couplings to quarks or leptons, is one of the simplest scenarios that 
can explain the dark matter. In this paper, we study in detail the impact of dark matter 
annihilation into the three-body final state WW (— > Wff') on the phenomenology of the 
inert doublet model. We find that this new annihilation mode dominates, in a relevant portion 
of the parameter space, over those into two-body final states considered in previous analysis. 
As a result, the computation of the relic density is modified and the viable regions of the model 
are displaced. After obtaining the genuine viable regions for different sets of parameters, we 
compute the direct detection cross section of inert higgs dark matter and find it to be up to 
two orders of magnitude smaller than what is obtained for two-body final states only. Other 
implications of these results, including the modification to the decay width of the higgs and 
to the indirect detection signatures of inert higgs dark matter, are also briefly considered. We 
demonstrate, therefore, that the annihilation into the three-body final state WW* can not 
be neglected, as it has a important impact on the entire phenomenology of the inert doublet 
model. 

1 Introduction 

Even though dark matter accounts for about 23% of the energy density of the Universe Q], we 
do not yet know what exactly it consists of. The identification of the dark matter particle is, 
indeed, one of the most challenging problems in astroparticle physics today. Over the years, many 
dark matter candidates have been proposed in different scenarios for physics beyond the standard 
model. Among them, the inert higgs -the lightest odd particle of the inert doublet model has 
earned a special place as a representative candidate of weakly interacting scalar dark matter. 

In the inert doublet model, a higgs doublet H2, odd under a new Z2 symmetry, is added to 
the standard model particle content. The lightest inert (odd) particle, H°, turns out to be stable 
and hence a suitable dark matter candidate. After being introduced in [2J, this model has been 
extensively studied in a number of recent works [3TU4|. 
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Recently, it was pointed out that the annihilation of dark matter particles can receive large 
contributions from three-body final states consisting of a real and a virtual massive particle [15HT7] . 
A fact that had been overlooked, notably, in all previous analysis of the inert doublet model. The 
annihilation into the three-body final state WW* (— > Wff), in particular, is expected to be 
important in the region M H o < M\y, which has already been shown to be viable and to feature a 
rich and interesting phenomenology [7]. In this paper, we revisit the inert doublet model of dark 
matter in view of these new contributions to the annihilation of dark matter particles. We will 
see that the inclusion of the three-body final state WW* not only modifies the viable parameter 
space of the model, but it also changes significantly the prospects for the detection of inert higgs 
dark matter. 

In the next section we present a brief introduction to the inert doublet model, followed, in 
section [3j by the calculation of the annihilation cross section of inert higgs dark matter into WW* 
(the analytical formulas are given in the Appendix). Section [4] is devoted to the analysis of the 
dependence of the three-body annihilation rate on the parameters of the inert doublet model. In 
section [5J the impact of the three-body final state on the inert higgs relic density is studied in 
detail. Then, we use those results, in section[Bl to derive the genuine viable parameter space of the 
inert doublet model. The implications of this new viable parameter space for the direct detection 
of inert higgs dark matter are investigated in section [7J Finally, in section [8j we concisely discuss 
other implications of the three-body final state regarding the indirect detection of dark matter 
and the decays of the higgs boson. 

2 The inert doublet model 

The inert doublet model is a simple extension of the standard model with one additional higgs 
doublet H 2 and an unbroken Z 2 symmetry under which H 2 is odd while all other fields are even. 
This discrete symmetry prevents the direct coupling of H 2 to fermions and, crucial for dark matter, 
guarantees the stability of the lightest inert particle. The scalar potential of this model is given 



where H\ is the standard model higgs doublet, and Xi and [ii are real parameters. Four new 
physical states are obtained in this model: two charged states, H ± 1 and two neutral ones, H° and 
^4°. Either of them could be dark matter. In the following we choose H° as the lightest inert 
particle, m 2 H0 < m^ ,m 2 H± , and, consequently, as the dark matter candidate. After electroweak 
symmetry breaking, the inert scalar masses take the following form 



by 



V =^il#i| 2 + l4\ H l\ + Ai|#i| 4 + A 2 |ff 2 | 4 + A 3 |i/i| 2 |i/ 2 | 2 
+ XilHlH.l 2 + ^ [{H\H 2 f + h.c] , 
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M2 + ^(A 3 + A 4 -A 5 )w 2 , 
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where v — 246 GeV is the vacuum expectation value of H\. Of pertinence to our study is the 
interaction term between a pair of H° and the higgs boson, which is proportional to Xl = (A3+A4 + 
As)/2. In addition to it, it is convenient to take m#o , Am^o = m^o — m^o , Am H ± — m H ± — m H a , 
and the higgs mass, m/,, as the remaining free parameters of the inert doublet model. 

In our analysis, we take into account all the known theoretical and experimental constraints 
on this model -see [I] and [7] . The requirement of vacuum stability imposes 

Ai,A 2 >0, A 3 ,A 3 + A 4 - |A 5 | > -2^\ 2 , (3) 

whereas the precise determination of the Z decay width at LEP requires that m^o + m#o > Mz 
and that m H + > Mz/2. Because of the specific decay modes considered in the analysis, the 
bound m H + > 79.3 GeV on the mass of a charged higgs from LEP [T8] can not be applied to 
this model. In [T5], the constraint m H + >70 — 90 GeV was derived using the results of the OPAL 
collaboration [20 . Some regions in the plane (rnjjo , m^o ) are also constrained by LEP II data, 
see [12]. Additionally, the inert doublet, H2, contributes to electroweak precision observables such 
as S and T. For the range of parameters we consider, however, compatibility with present data 
is easily achieved. Finally, from section [5] on, we require also that the relic density of inert higgs 
dark matter be compatible with the observed dark matter density. 

In previous works [4l[51[7 HT^l[2Tj several aspects of this model, including constraints from 
present data and prospects for dark matter detection, were studied. It turns out that the dark 
matter constraint can only be satisfied for restricted values of mjyo . Three viable regions can be 
distinguished: a small mass regime with m H o ~ 8 GeV [6"1[TT] . a large mass regime with m# > 500 
GeV [7,22,23 and an intermediate mass regime (mjj Q < Mw) [HIT] . This intermediate mass regime 
gives rise to a very rich phenomenology [7_ , with significant direct detection cross section and good 
indirect detection prospects via gamma rays [5]. It also coincides with the mass range where the 
annihilation into the three-body final state WW* , a process not considered in any previous work 
on this model, is expected to be particularly relevant. In this paper, we take into account, for the 
first time, the impact of the WW* final state on the phenomenology of the inert doublet model. 

3 H°H° annihilation into WW* 

It has been recently emphasized that dark matter annihilations could receive large additional 
contributions from three-body final states consisting of a real and a virtual massive particle [T5] . 
The inert higgs dark matter, in particular, could annihilate into WW*, ZZ* , hh* , and ti* . The 
latter two, however, are certainly irrelevant. The branching ratio into them is small and, in any 
case, the region where they could give a significant contribution (for m#o just below the h and t 
thresholds) is not consistent with the relic density constraint. In that region, the annihilation into 
W + W~ and ZZ is so effective that it drives the thermal abundance of inert higgs dark matter well 
below the WMAP measurement. Including the additional annihilation into hh* and ti* would, 
at most, yield an even smaller relic density. That leaves us only with the gauge bosons three- 
body final states. They seem promising because the region below the W and Z thresholds partially 
overlaps with the intermediate viable region of the inert doublet model, M H a < Myy. Between the 
two, we expect the WW* final state to be more important, for we know that its real counterpart 
dominates the annihilation branching ratio in the region m#o > My/,Mz- Besides, the virtual 
Z is farther off-shell than the virtual W within the viable region of the model, m#o < Mw, 
giving an additional suppression. It turns out that the ZZ* final state is actually negligible with 
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Figure 1: The Feynman diagrams that contribute, in the unitary gauge, to the H H annihilation 
into the three-body final state WW* — > Wff within the inert doublet model. For the H + - 
mediated diagram the exchange diagram (not shown) must also be taken into account. 



respect to WW* over the entire viable parameter space of the inert doublet model. Thus, the only 
three-body final state that can modify the phenomenology of the inert double model is WW* . 

Three different diagrams contribute to the annihilation of dark matter into WW* (H°H° —> 
WW* — > Wff) in the inert doublet model -see figure [TJ The amplitude of the first diagram 
depends only on gauge parameters whereas the second and the third also depend respectively on 
rriH+ and on Az, and rrih- The contribution from the direct annihilation diagram (left diagram) 
is usually dominant whereas the one from the iJ^-mediatcd diagram is typically small. The 
higgs-mediated contribution can be important, particularly close to the higgs resonance, where 
it becomes dominant. a(H°H a — > WW*) is of utmost importance in our discussion, as it enters 
explicitly into the indirect detection rates and into the computation of the relic density. The 
analytical result for this cross section can be found in the Appendix [A] (t(H°H° — » WW*) 
depends very weakly on m^o (only through the higgs width), and on m#+ (the H + mediated 
diagram is suppressed by the t-(u-)channel propagator). So, we will study its dependence on 
m H o, Al (sign and magnitude), and rrih- 

The two-body annihilation rate, on the other hand, is determined by higgs-mediated processes 
into light fermions. So, it is proportional to \ 2 L and it is dominated by the bb final state. In 
spite of being formally of higher order, the three-body process can compete with the two-body 
ones thanks to the Yukawa suppression present in the latter and to the large multiplicity of final 
states associated with WW* (— > J2f Wff). To check that it is indeed the case, we must compute 
the three-body annihilation cross section, a(H°H a —> WW* — > Wff), and compare it with the 
two-body one. 

To begin with, let us study the behaviour of the annihilation rate at low velocities, crv, with 
respect to the parameters of the model. Figure [5] compares the two-body and the three-body 
(WW*) annihilation rate as a function of m H o for two different higgs masses, 120 GeV (left 
panel) and 150 GeV (right panel), and the two possible signs of A^. From the left panel, we see 
that, as expected, the three-body cross section generically increases as m#o gets closer to Mw 
The atypical behaviour observed around the higgs resonance, m#o ~ 60 GeV, is the result of the 
interference between the purely gauge diagram and the higgs mediated diagram, as explained in 
Appendix iBl Because of such interference, the three-body cross section for A^ > (dash-dotted 
line) is larger than that for Al < (dashed line) above the higgs resonance but smaller than it 
below the resonance. In any case, the crucial point for us is that the three-body cross section 
is not negligible at all. It becomes larger than the two-body one for mno > 62 GeV or 67 GeV 
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Figure 2: Comparison between the three-body and the two-body annihilation rate, crv, as a 
function of the dark matter mass for the two possible signs of A^. In the left panel rrih = 120 
GeV whereas in the right panel m,h = 150 GeV. The other parameters were taken as Amj = 
Am H ± = 50 GeV, |A L | = 10" 2 . 



depending on the sign of A^. 

For mil = 150 GeV (figure[2j right panel), the effect is even more pronounced. In this case, the 
three-body cross section dominates the annihilation rate in almost the whole range m^o > 50 GeV 
-a fact partially due to the suppression of the two-body annihilation that is expected for a higher 
higgs mass. We see that the only regions where the two-body cross section is larger are two 
narrow mass intervals around the resonance where the interference effects between the higgs and 
pure gauge contribution suppress the three-body cross section. Also notice that right at the 
resonance, the two-body annihilation rate is larger than the three-body one for to/, = 120 GeV 
whereas it is the other way around for rrih = 150 GeV. This is in agreement with the known 
result that in the Standard Model a 120 GeV higgs boson decays dominantly into two-body final 
states whereas a 150 GeV higgs boson decays mainly into three-body final states. For even higher 
values of the higgs mass, the differences between the Xl > and Xl < cases will tend to fade 
out as the higgs resonance moves further away from the relevant parameter space. We have thus 
illustrated, via figure [2j the importance of cr(H a H a — > W / M / *)v and its dependence on m^o, rrih 
and the sign of A^ . 



4 Effects on av - a systematic analysis 

With the goal of understanding the relevance of the three-body final state, in this section we 
perform a systematic analysis of the ratio between the three-body cross section a(H°H° -> WW*) 
(sometimes denoted simply as a(WW*)) and the two-body one a(H°H° — » Ylf ff) (denoted also 
as cr(2-body)) within the inert doublet model. For the higgs mass, we consider three typical values 
compatible with electroweak precision data: 120, 150, and 200 GeV. 

Figure [3] displays the ratio of cross sections as a function of m#o for different values of the 
higgs mass. For this figure Xl was set to 10 -2 whereas Am^o and Am^+ , which hardly affect the 
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results, were set to 50 GeV. From the figure we see that, the ratio is larger for m/j = 200 GeV over 
the whole range of m#o. For that higgs mass, the ratio increases with m H o and it is larger than 
1 in the entire range My/ > m H o > 50 GeV, implying a dominance of the three-body final state 
over the two-body ones. For the other two higgs masses, we see that the ratio tends to increase 
with m H o but it features a narrow dip before reaching the higgs resonance. For = 120 GeV the 
three-body final state dominates the cross section for mjyo > 62 GeV whereas for rrih — 150 GeV 
it dominates it over the range to#o > 52 GeV except for a small mass range around 73 GeV. 
Independently of the higgs mass, the ratio can be larger than 100 close to the W threshold, 
indicating a strong dominance of the three-body final states in that region. 

The annihilation cross section into WW* is affected by Al via the higgs mediated diagram. 
This dependence is illustrated in figure 21 which shows as a function of m#o , the ratio of cross 
sections for two different values of Al (10~ 2 and 10~ 3 ) and two different higgs masses (120 GeV 
and 150 GeV). As expected, a smaller Al suppresses the two-body cross section, which being 
higgs mediated is proportional to X L , increasing the relevance of the three-body final state. In 
fact, the ratio is about two orders of magnitude larger for Al = 10~ 3 than for Al = 10~ 2 . As 
a result, for Al = 10 -3 the three-body cross section is larger than the two body one essentially 
over the entire m#o range we consider, and the ratio can reach values above 10 4 close to the W 
threshold. Another feature observed in the figure is the displacement of the position of the dip, 
which due to the smaller interference effects moves closer to the resonance for smaller values of 
Al -see appendix 151 

As already shown in figure [2j the three-body cross section depends also on the sign of Al, due 
to the interference between the higgs mediated amplitude, which goes like Al, and the other two 
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Figure 4: This figure illustrates the dependence on of the ratio between the three-body and the 
two-body annihilation rate. It shows cr(3-body)/cr(2-body) as a function of mjyo for two values of 
Xl (10 -2 and 10 -3 ) and two different higgs masses (120 and 150 GeV). For the other parameters, 
we take Am^ = Amjy± = 50 GeV. 



contributions. Figure [5] shows the ratio of cross sections as a function of m#o for the two possible 
signs of Xl and two different higgs masses: 120, 150 GeV. The other parameters were taken as 
|Al| = 10~ 2 , Am^o = Am H + = 50 GeV. The sign of Xl clearly influences the three-body cross 
section, particularly around the higgs resonance, where the interference is stronger. Notice that the 
dip moves across the resonance when Xl changes sign. Moreover, the curves for Xl > and < 
cross each other exactly at the resonance, as expected. Indeed at that point, the annihilation cross 
section is entirely determined by the higgs mediated diagram, so cr(WW*) cx X\ there. Notice 
that the ratio is larger for Xl > than for Xl < above the resonance, whereas it is the other way 
around below the resonance. Thus, the largest effect will be obtained for Xl > if m/, = 120 GeV 
but for X L < if m h = 150 GeV. 

This detailed study of a(H°H° WW*)/a(H°H° //') clearly demonstrates that the 
relevance of the annihilation into the three-body final state WW* is a generic feature of the inert 
doublet model. No fine-tuning is necessary to find large effects. In the next section we investigate 
the implications of this new process on the calculation of the inert higgs relic density. 

5 Effects on the relic density 

To compute the dark matter relic density, one must solve the following Boltzmann equation, 

% = \f^pM p (av)(Y(Tf - Y eq (Tf) , (4) 
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Figure 5: This figure illustrates the dependence on the sign of Xl of the ratio between the three- 
body and the two-body annihilation rate. It shows er(3-body)/cj(2-body) as a function of m^o 
for the possible signs of Xl and two different higgs masses (120 and 150 GeV). For the other 
parameters, we take Am^o = Amj± = 50 GeV and |Al| = 10~ 2 . 



where Y(T) is the dark matter abundance, defined as the number density divided by the entropy 
density, g* is an effective number of degrees of freedom, M p is the Planck mass, and Y eq (T) is the 
equilibrium thermal abundance, (crv) is the thermally averaged annihilation cross section, which 
must include all relevant annihilation and coannihilation processes. It corresponds to 



where i,j run over all annihilating and coannihilating particles, gi is the number of degrees of 
freedom of particle i, rrii is its mass, pij is the momentum of the incoming particles in the center 
of mass frame, and aij-^i is the total annihilation cross section of particles i, j into Standard Model 
particles k, I. Integrating equation ((4]) down to today's temperature T = To leads to the present 
dark matter abundance, Y(T ). From it, the relic density can be obtained as 



where M^ m is the dark matter particle mass. This standard approach for the calculation of the 
relic density has been implemented in publicly available software such as DarkSUSY 24, 25 and 
micrOMEGAs [2BH25]. 

To properly compute the inert higgs relic density we need to include in equation ([5]), in addition 
to the usual two-body annihilation and coannihilation processes, the annihilation into the three- 
body final state WW*. That is the only modification we need to make to the above procedure 
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Figure 6: Comparison between the 3-body and the 2-body relic density, f2/i 2 , as a function of the 
dark matter mass for the two possible signs of Al. In the left panel rrih = 120 GeV whereas in 
the right panel rrih = 150 GeV. The other parameters were taken as Am^o = Am H ± — 50 GeV, 
|A L | = 10- 2 . 



-equations ((4]) and ([6]) remain the same. For our calculations, we have used a modified version of 
micrOMEGAs, in which we incorporated the annihilation into the three-body final state WW* in 
the evaluation of ([5]). That way, we can accurately compute the relic density of inert higgs dark 
matter including 3-body final states. 

In the following we compare the relic density obtained for two-body final states only (denoted 
as SI (2-body)) with that predicted including also the final state WW* (denoted as (3-body) and 
referred to as the 3-body relic density) for different values of the parameters of the inert doublet 
model. 

Figure |6] shows the relic density derived including only 2-body processes (denoted 2-body) and 
including 2-body and 3-body processes relic density (denoted 3-body) as a function of m H o for the 
two possible signs of Al and two different higgs masses: 120 GeV (left panel) and 150 GeV (right 
panel). The remaining parameters were taken as |Al| = 10 -2 , Am^o = Am H + = 50 GeV. First 
of all, notice that the two-body relic density (solid and dotted lines) depends weakly on the sign 
of Al, and only for m^o close to Myy- This behavior is to be expected as it is the annihilation 
into W + W~ (both being real particles) that brings such a dependence into play, and it is only 
for m#o close to M\y that such annihilation can take place in the early Universe. The 3-body 
relic density is observed to be always equal or smaller than the two-body one -as anticipated- 
and to depend on the sign of Al- For rrih = 120 GeV (left panel) the effect of the 3-body relic 
density is significant above the resonance, and the predicted relic density is smaller for Al > 
(dash-dotted line) than for Al < (dashed line). For rnj^a close to M\y, the dark matter particles 
may have enough kinetic energy in the early Universe to annihilate into W + W~; consequently, 
the effect of the three-body final state becomes less relevant in that region. For mn — 150 GeV 
(right panel), the 3-body relic density is significantly smaller than the two-body one over the entire 
mass range we consider. Regarding the sign of Al, both signs give approximately the same relic 
density below the resonance and differ from each other above it. These two figures demonstrate 
that the three-body final state WW* affects in a relevant way the predicted relic density of inert 
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Figure 7: Figure illustrating the dependence in the higgs mass of the ratio between the relic density 
including the three-body final state and the relic density for two-body final states as a function 
of the dark matter mass. The other parameters were taken as Ato^o = Am H ± = 50 GeV and 
A L = 10- 2 . 



higgs dark matter. 

To better illustrate the effect of the three-body final state, we next study the ratio between 
the three-body and the two-body relic density for different sets of parameters. Figure [7J shows 
this ratio as a function of mjjo for different values of the higgs mass. A ratio equal to 1 means 
that the three-body process gives a negligible correction to the calculation of the relic density. 
Clearly, that is not the case. The ratio tends to 1 for m H o close to M\y, where the annihilation 
into W + W~ is efficient, and for m H o <C Mw, where the three-body annihilation is suppressed, 
but in the intermediate region the three-body final state plays a major role, giving rise to a correct 
relic density significantly smaller than the two-body one. An effect that, as observed in the figure, 
is present independently of the higgs mass -although its precise magnitude will certainly depend 
on Trifi . Notice, from figure [71 that the two-body approximation may overestimate the predicted 
relic density by more than one order of magnitude. Moreover, a significant deviation from the 
two-body result can take place over a wide range of rn//o; depending on the higgs mass, it could 
extend more than 30 GeV below M\y. 

The mass of the CP-odd scalar, m^o, may affect the H° relic density via coannihilation pro- 
cesses. If the mass splitting between H° and A , Am^o, is small, the process H°A° — >• Z* — > //' 
increases the annihilation rate and help reduce the relic density. Let us point out that, in the 
inert doublet model, there is no need to consider possible coannihilations into three-body final 
states because the two-body ones, being gauge processes, are unsuppressed. Since coannihilations 
increase the annihilation rate, they are expected to reduce the importance of the three-body final 
state. That is exactly what is seen in Figure which shows the effect of a smaller Am^o on the 
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Figure 8: Figure illustrating the dependence on Amjo of the ratio between the relic density 
including the three-body final state and the relic density for two-body final states only. Arra^o 
affects £lh 2 only through coannihilations effects, which are important for small mass splittings 
(Am A a — 10 GeV) but not for large ones (Am^o = 50GeV). Two different values of the higgs 
mass are shown. The other parameters were taken as Am H ± = 50 GeV and Xl = 10~ 2 . 



relic density. It displays the ratio between the 3-body and the 2-body relic density for two values 
of Am^o and two different higgs masses. When Am^o = 10 GeV, coannihilations are important 
(dotted- and dashed-lines) , and the three-body final state is less relevant -the ratio is closer to 1— 
than for Am^o = 50 GeV when coannihilations are suppressed (solid- and dashed-dotted lines). 
For to /j = 150 GeV the deviation due to the different values of Am^o is significant for m H o below 
65 GeV whereas for mh = 120 GeV it is so for m#o between 60 GeV and 75 GeV. Notice also 
that close to the higgs resonance, the H°H° annihilation rate is enhanced, making coannihilation 
effects negligible. That is why the curves for Am^o = 10 GeV and Am^o = 50 GeV coincide in 
the region m#o < Although coannihilation effects slightly reduce its relevance, the effect of 

the three-body final state remains significant over the parameter space of the inert doublet model. 



6 The genuine viable parameter space 

Instead of computing the dark matter density for a given set of parameters, we are usually inter- 
ested in using the relic density as a constraint on the parameter space of the model. The viable 
parameter space is determined by requiring that the predicted relic abundance be compatible 
with the observed density of dark matter [T] . For the inert doublet model, this viable parameter 
space had been obtained in previous works, but using the two-body relic density, which, as shown 
in the previous section, is not a good approximation in the intermediate mass regime. In this 
section, the genuine viable parameter space of the inert doublet model, obtained by including the 
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Figure 9: The viable parameter space for mh — 120GeV. Along the lines ilh 2 = 0.11. The thick 
lines are the result including the final state WW* , the thin lines correspond to 2-body final states 
only. 



three-body annihilation into the computation of the relic density, is derived. We study in detail its 
dependence on the parameters of the model and we demonstrate that it is significantly different 
from that found for two-body annihilations. 

For dcfinitcncss, wc focus on the following three interesting cases: = 120 GeV with Az, > 0, 
rrih — 150 GeV with Az, < 0, and mt = 200 GeV with Az, < 0. Figure [9] shows the viable parameter 
space of the intermediate mass range of the inert dark matter model in the plane (Az,,mzzo) for 
rah — 120 GeV, Am H ± — 50 GeV, and two different values of AmA , 10 GeV and 50 GeV. These 
two values of Am,A are chosen so as to indicate the possible effect of having (Amyi = 10 GeV) or 
not (AmA = 50 GeV) significant coannihilation processes. The thin lines in this figure correspond 
to the viable regions if only two-body final states are considered. The thick lines, on the contrary, 
correspond to the genuine viable regions, those obtained by taking into account two- and three- 
body final states in the calculation of the relic density. Two important results are clearly observed 
in this figure. First, the value of Az for a given mass may be substantially smaller once three- 
body final states are taken into account. The difference between the two Az, associated to a fixed 
value of Arn^o and m H a could amount to one order of magnitude. Second, the viable parameter 
space shrinks toward lower masses. When only two-body final states are considered, the maximum 
values of mzzo allowed are about 75 and 76 GeV respectively for Am^o = 10 and 50 GeV. Once the 
three-body final state is included, the maximum mzzo moves respectively to about 71 GeV and 73 
GeV. This reduction of the viable mass range is entirely due to the effectiveness of the three-body 
annihilation, which drives the relic density below the WMAP bound for mzzo > 71, 73 GeV. In 
that region, no value of Az, is allowed, for it is the direct annihilation (figure [U left diagram), a 
gauge process independent of Az, that determines the relic density. As observed in this figure, the 
difference between the two-body and the three-body viable parameter space is quite significant, 
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Figure 10: The viable parameter space for mh = 150GeV (left panel) and mh — 200GeV (right 
panel). Along the lines ilh 2 = 0.11. The thick lines are the result including the final state WW*, 
the thin lines correspond to 2-body final states only. 



independently of possible coannihilation effects. 

Similar results are found for other higgs masses, as illustrated in figure [TU1 for m/j = 150 GeV 
(left panel) and = 200 GeV (right panel). In both cases < was considered and the 
conventions used were the same as in figure 1 As a consequence of the three-body final state 
contribution to the annihilation rate of inert higgs dark matter, the required value of Xl is smaller 
at any given mass, and the maximum allowed value of mjjo gets reduced by several GeVs. 

Figures [9l and [T0l confirm that the modification of the viable parameter space, induced by the 
annihilation into the three-body final state WW*, is a generic feature of the inert doublet model. 
A feature that is present over a wide range of mjj° independently of the other parameters of the 
model. It is precisely because of this generality that the intermediate mass regime of the inert 
doublet model must be revisited, as we do in this paper, in view of these new processes. 

These modifications to the viable parameters space are of crucial relevance because they affect 
other phenomenological aspects of the model. In fact, the first step in the analysis of a given dark 
matter model is usually the determination of its viable parameter space. Once it is obtained, one 
can study the specific signatures or predictions of the model within such viable regions. If these 
are changed so are its predictions and signatures. In the next section we show that, as a result 
of the new viable parameter space, the direct detection cross section of inert higgs dark matter is 
considerably reduced. 

7 Direct detection of inert higgs dark matter 

Dark matter can be detected via elastic scattering with terrestrial detectors, the so-called direct 
detection method. From a particle physics point of view, the quantity that determines the direct 
detection rate of a dark matter particle is the dark matter-nucleon scattering cross section. In 
the inert higgs model, the H°N scattering process relevant for direct detection is higgs-mediated, 
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Figure 11: The direct detection cross section along the viable regions of the inert doublet model 
for nih = 120 GeV and Am fl + = 50 GeV and two different values of Am^o. 



with a cross section, <Jh n, given by 



where m r is the reduced mass of the system, vtln is the nucleon mass that we took equal to the 
proton mass, and / is the nucleon form factor, taken equal to 0.3 for the subsequent analysis (see 
e.g. the discussion in QT]). Hence, (Jh°n is proportional to \ 2 L . Given the new allowed values of 
\l that were derived in the previous section, we foresee that a H o n will be significantly reduced 
with respect to the two-body result used, until now, in the literature. 

Figure [TT] shows the prediction for <Jh°n along the viable lines of the inert doublet model for 
rrih — 120 GeV. As before, thin lines correspond to the two-body result while thick lines to the 
three-body one. For comparison, in this figure we also show, as a dotted line, the current limit 
from CDMS [30 . Notice from the figure that the correct direct detection cross section can be 
more than two orders of magnitude smaller than the one obtained for two-body final states. In 
particular, for rrih — 120 GeV the genuine direct detection cross section turns out to be well below 
present bounds. 

Analogous results follow also for other values of the higgs mass. Figure [T2l shows the same cross 
section but for rrih = 150 GeV (left panel) and rrih = 200 GeV (right panel). In both we observe 
that, over a wide range of mfjo, the correct cross section (denoted by 3-body) is significantly 
smaller than the one derived after including only 2-body processes in the computation of the relic 
abundance (denoted by 2-body). Notice, for instance, that for rrih — 200 GeV the region above 
the CDMS bound moves from m#o < 71 GeV to m^o < 62 GeV when the three-body final state 
is taken into account. 




(7) 
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Figure 12: The direct detection cross section along the viable regions of the inert doublet model 
for nih — 150 GeV (left panel) and mt — 200 GeV (right panel). For the other parameters, we 
took Am H + = 50 GeV and Am A o = 10, 50 GeV. 



The inclusion of three-body final states is, therefore, mandatory if one wants to make mean- 
ingful predictions on the prospects for the direct detection of inert higgs dark matter. 

8 Other implications 

In this section we briefly address other possible implications of the three-body process on the phe- 
nomenology of the inert doublet model. Specifically, we show, in section [8TTT that the annihilation 
branching ratios can be dominated by the three-body final state, with important implications for 
the indirect detection of dark matter, and we discuss, in section [8?2| the modifications to the decay 
width of the higgs boson. 

8.1 Indirect detection 

The indirect detection signals of inert higgs dark matter are also altered by the existence of the 
three-body final state WW* . On the one hand, these signals should be now computed along new 
regions, due to the modified viable parameter space. On the other hand, in these new regions 
the annihilation cross section and branching ratios typically receive large corrections from the 
three-body final state WW* . Indeed, we already saw, in section 21 that the annihilation of inert 
higgs dark matter may be dominated by three-body final states, rather than by the two-body final 
states considered in previous works. As a result, the spectrum of photons, neutrinos, positrons and 
antiprotons expected from inert higgs annihilation will be different, changing its indirect detection 
prospects. Even the gamma ray lines from the one-loop annihilation into two photons [8] will be 
affected by these new contributions. A detailed study of the implications of three-body final 
states for the indirect detection of inert higgs dark matter is beyond the scope of the present 
paper. Here, we just want to demonstrate that, contrary to what has been assumed in earlier 
analysis on indirect detection of inert higgs dark matter [10], bb is not necessarily the dominant 
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annihilation channel for m H a<Mw ■ The three-body final state WW* turns out to be dominant 
over a sizeable region of the viable parameter space. 

Figure [13] shows the annihilation branching ratio into the three-body final state WW* as a 
function of m#o along the viable regions of the inert doublet model. Each line corresponds to a 
given value of ran and a given sign of \l- For this figure we set Am^o = Ato#+ = 50 GeV but 
the results are similar for other allowed values. Notice that the branching into WW* is indeed 
significant: it amounts to more than 10% for m H o > 55 GeV independently of the higgs mass. 
Moreover, it reaches values close to 1 for 55 GeV < m H a < 60 GeV if rrih — 120 GeV and also for 
65 GeV < rriffo < 72 GeV independently of the higgs mass. In view of these results, the indirect 
detection signatures of inert higgs dark matter will have to be revised. In a future work, we plan 
to carry out such an analysis. 

8.2 Higgs decays 

In the inert doublet model, the higgs boson can decay also into H°H° and A A , increasing the 
higgs decay width and modifying its branching ratios. A result that is of great interest for higgs 
searches at colliders [9]. Here we simply illustrate how the higgs decay width is modified when 
including the three-body final state in the determination of the relic abundance. As before, we 
assume that H° accounts for the dark matter of the Universe. 

The contribution to the higgs decay with from the decay into the inert scalars is proportional 
to A|, so it will be affected by the three-body final state WW* via the new viable parameter 
space. Figure [14] shows, as a function of rn^-o, the ratio between the higgs decay width in the inert 
doublet model and that in the Standard Model, T(h)/T sm (h), for = 150 GeV. Obviously, 
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Figure 14: Ratio between the higgs branching ratios in the inert doublet model and in the standard 
model along the viable regions for mt = 150 GeV. Xl was taken to be negative and Am H ± = 50 
GeV. 



the higgs width is only affected for m#o below rrih/2- For the inert doublet model we consider 
the result with and without taking into account the three-body final state in the calculation of 
the relic density, and we set Am#+ = 50 GeV, < 0. First of all, notice that the higgs decay 
width can indeed be much larger than what is predicted by the Standard Model. Second, the 
higgs width in the IDM is slightly smaller once three-body final states are taken into account -see, 
for instance, the difference between the solid and dotted line or between the dash-dotted and the 
dashed line. Notice, however, that the region where the deviation from the standard model result 
is significant, for 50 GeV < m H o < 60 GeV, is partially excluded by the present direct detection 
bound from CDMS -see figure [TJ] This result emphasizes the importance of doing a consistent 
analysis of the inert doublet model, one that simultaneously takes into account all the relevant 
effects and the different constraints that can be imposed on the model. 

9 Conclusions 

We studied the impact, on the phenomenology of the inert doublet model, of dark matter annihi- 
lation into the three-body final state WW* . After analyzing its dependence on the parameters of 
the model, the annihilation cross section into WW*, a(H°H° — > WW*), was shown to dominate 
the total dark matter annihilation cross section over a relevant portion of the parameter space. In 
consequence, the predicted relic density differs considerably from that found in earlier works. We 
examined in detail the dependence of the inert higgs relic density on the higgs mass and on Al, 

Notice that T(h) and T S m(h) includes the h — > WW* channel in all the curves 
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as well as the possible role of coannihilations. The genuine viable parameter space of the inert 
doublet model was derived for different values of the higgs mass and of Am^o , emphasizing the 
differences with respect to the two-body results. We also investigated the direct detection cross 
section of inert higgs dark matter and found that it can be more than two orders of magnitude 
smaller than what is predicted without including the three-body final state. Finally, we briefly 
consider some implications of these new annihilation processes on the decay width of the higgs 
boson and on the indirect detection of inert higgs dark matter. Summarizing, the inclusion of the 
three-body annihilation is mandatory, as it strongly affects the entire phenomenology of the inert 
doublet model. 
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A Analytical formula for a(H°H° — > WW*) in the inert dou- 
blet model 

Here we provide the total amplitude squared for the process H°H° —> WW*, which has been 
directly obtained with the Calchep package [31] . The relevant Feynman diagrams are shown 
in figure [1] We denote by A4 P the amplitude for the direct annihilation diagram (left diagram 
in figure [T]) and by M. s the higgs mediated diagram (right diagram). M u .t correspond to the 
H + mediated diagrams (middle). In the following, p\ and p 2 denotes the 4- momentum of the 
annihilating Hq, p 3 and p 4 are the momentum of the 2 fermions produced in the decay of the 
virtual W* (pw* = P3 + Pi) and p$ is the momentum of the real W . 

, M , 2 _ 9 6 - 2m H + + 2(-p 1 .p 2 +Pl.p 3 +Pl.p4+P2-P3+P2-P4) 

' Pl ' 4M^, ((( p 3+ p 4 )2_ M 2 /) 2 + (rH/MH/) 2 ) 

-o 6 

\Mt\ 2 = &M 2 Df KKfo - ™H+? - A&(M^ + 4(-pi.p2 +P2-P3 +P2-P4))) 

(2m% + 4p 1 .p 3 p 1 .p 4 ~ m 2 Ho (Mw + 2(-p 1 .p 2 + Pi-P3 + Pi-Pa +P2-P3 + P2-Pa)))] 

D t = (t-m 2 H+ ) 2 ^(((p 3 + Pi ) 2 -M^) 2 + (r w M w ) 2 ) 

|A/,tl ' 2 = 8M 2 £>J ~ ((m "° ~ m ^ +)2 ~ M w( M w+<-Pi-P*+Pi-P3 +P1-P4))) 

(2m% + 4:P2-P3P2-P4 - m 2 H { M w + 2 (-Pi-P2 + Pi-P3 + Pi-Pa + P2-P3 + P2-Pa)))] 

D u = {u-m^Y^iiipz+pif-MlrY + iTwMw) 2 ) 
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g 6 

|A^ S | 2 = 2 [(ffij, - M2) 2 ( 2 ^H + + 2(pLP2 - P1.P3 - P2-P3)(P1-P2 - Pi-Pi - P2-Pi) 

* V *W P 

+M^r(-3p 1 .p 2 + 2{pi.p 3 +Pl.p 4 +P2-P3 +P2-P4)) 
-m 2 Ha ('AMw + 2(-2p 1 .p 2 +pi.p3 + P1-P4 + P2-P3 + P2-Pi)))\ 

D s = ((s-ml) 2 + {m h T h ) 2 )*{{{p3+ Pi ) 2 - M 2 ,) 2 + (T W M W ) 2 )) 

2M s Ml = [K -m)(m Ho +M2) 

lvl w P s 

(2m 2 Ho - - 2(-p 1 .p 2 + P1-P3 + P1-P4 + P2-P3 + P2-P4))] 

D sp = (s-m 2 h )*(((p 3 +P4) 2 - M 2 ,) 2 + (T W M W ) 2 ) 
2MsM * = 2M^ st [(m ^~^ )(2m6g « 

+™ffo(~ 2m ff+ + M VV + 4pi-P2 - 2(pi.p 3 +P2-P3 +P2-P4)) 

+m^+(M^pi.p 2 - 2(pi.Pa +P1-P4P2-P3 

+P1.P3(2P1.P4 +P2-P4) ~Pl-P2(Pl-P3 +P1-P4 +P2-P3 +P2-P4))) 

+M^(M^pi.p2 + 2(-pi.j^ + 2pi.p3Pi.pi 

~Pl-PiP2-P3 ~P\-P3P2-P\ +Pl-P2(Pl-P3 + Pl-PA +P2-P3 +P2-Pi))) 
+ m H (- M W - Mw(Pl-P2 + 2(P!.P3 + P!.p 4 + P2-P3 + P2-P4)) 
+2(pi.pa +P1-P4P2-P3 + Pl-P3(2pi-P4 +P2-P4) 
-pi-P2{Pl-P3 +P1-P4 +P2-P3 +P2-P4)) 

+m 2 H+ (Mw + 2(-2p 1 .p 2 +P1-P3 +P1-P4 +P2-P3 +P2-P4,))))] 

D s t = (a - m£) * (((pa + p 4 ) 2 - M^) 2 + (%M ff ) 2 ) * (t - m 2 H+ ) 
2MsM " = 2M^Z? SM [ ~ (m ^~^ )(2m ^ 

+ m Ho( _2TO H+ + M VV + 4 PlvP2 - 2(pi.p 3 +Pl.p 4 +P2-P3 +P2-P4)) 

+m 2 H+ (MwPi-P2 - 2(pi.pl +P1-P4P2-P3 

+ (Pl-P3 + 2p 2 .p3)P2-P4 ~Pl-P2(Pl-P3 +P1-P4 + P2-P3 + P2-Pi))) 
+Mw(MwPl-P2 + 2{-p\.p\ ~P\.P4P2-P3 

-Pi -P3P2 -Pi + 2p 2 -P3P2 -Pi+Pl -P2 (Pl-P3+Pl -Pi + P2 -P3 + P2-Pi))) 
+ m H„(- M W - Mw(Pl-P2 + 2(pi-P3 + P1-P4 + P2-P3 + P2-P4)) 
+2(p 1 .p\ +P1-P4P2-P3 + (Pl-f>3 + 2p 2 .p3)P2-P4 
~Pl-P2{P\-P3 +P1-P4 +P2-P3 + P2P4)) 

+m 2 H+ (Mw + 2(-2p 1 .p 2 +P1-P3 +P1-P4 +P2-P3 +P2-P4,))))] 

D su = (s-m 2 )*(((p3+P4) 2 -M^) 2 + (T w M w ) 2 )*(u-m 2 H+ ) 

a 6 

2M P M\ = — | - [-2m% - M^ Pl .p 2 + m 2 Ho (M^ + 2( Pl .p 3 + Pi.p 4 + P2-P3 + P2-Pa)) 

iM w D pt 

+ 2(p 1 .pj +Pi.p 4 P2-P3 + Pl-P3(-2p 1 .p 4 + p 2 .p 4 ) 
~Pl-P2{Pl-P3 + P1-P4 + P2-P3 + P2-P4))} 
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D P t = (((P3+P4) 2 - M 2 ,) 2 + {T w M w ) 2 )*(t-m 2 H+ ) 
2M p Mi = 2 [2m% + Mlr Pl .p 2 - m 2 Ho (M^ + 2( Pl .p 3 + Pl . Pi + p 2 .p 3 + P2 . Pi )) 



+ 2(-pi.p2 - P1-P4P2-P3 - P1-P3.P2-PA + 2 P2 - P 3P2-P4 
+Pl-P2{Pl-P3 +P1-P4 +P2-P3 +P2-PA))] 

D P u = {{{ P 3+Pi) 2 ~ M 2 ,) 2 + {T w M w f)*{u-m 2 H+ ) 

^6 



2M * M l = . . J n [- (m^ + m% + - 2mf Io (m 2 H+ - 2M^) 



t - -g 

-Mw(Mw + 2(p!.p 3 +pi.p 4 +P2-P3 + P2-Pi))){2m 2 Ha pi.p 2 - M^ v p 1 .p 2 

+ 2{px.pl +P1.P4P2-P3 + P1-P3P2-P4 - Pl-P2(Pl-P3 + P1-P4 + P2-P3 + P2-P4,)))] 



D tu = (t - jti 2 h+ ) * (((p 3 + p±) 2 - M^f + {T w M w ) 2 )*{u-m 2 H+ ) 

The cross section a(H°H° —> WW*) can then be obtained, from the total amplitude, in the usual 
way. 

B Interference effects around the higgs resonance 

Around the higgs resonance, cancellations between different contributions clearly appear in cri>(3-body). 
They take place at mn < mh/2 for \l > and at rriH > mh/2 for Al < 0. This fact can be 
easily understood by looking at the square of the sum of the two relevant amplitudes, which are 
the one of the point like diagram Ai p and the one of the higgs mediated diagrams M s . They can 
be written as 



Xlv 2 9 2 
ml - imhTh 



iM s = £—2 —O w »0$y. (8) 



where, in a similar way to |16) , we have defined Own as the wave function of the outgoing real W 
and Ow n as the contribution from the virtual gauge boson and the two outgoing fermions: 

= t : ^ 772 — — j=u 3l v±. (9 

(P3 +P4) 2 - M w ~ ilwMw y/2 



For the square of their sum we obtain 

\V«p + * M *\ 2 - t (7^g^4l [(s " ml) + 2Xlv2? (10) 



Since s ~ 4m^ o at low velocities, a negative interference between the two terms should occur 
for mji < mh/2 if Xl > and for mu > mh/2 if A^ < 0. That explains why the 3-body 
annihilation cross section is larger for Al < than for Al > if m^a < mh/2 but it is the other 
way around for m^o > mh/2. Moreover, a cancellation between these two contributions takes 
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place at \m,H — nih/2\ ~ 1/2\Xl\v 2 /nih- For = 120,150 GeV and Xl = 10~ 2 , that position 
correspond to |m# — m,h/2\ ~ 2.5,2.0 GeV, in good agreement with figures [3j [4] and [5j which 
were obtained integrating numerically the cross section including all the contributions displayed 
in appendix IA1 
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